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We study high-density electron-hole (e-/i) systems with the electron density slightly higher than 
the hole density. We find a new superconducting phase, in which the excess electrons form Cooper 
pairs moving in an e-h BCS phase. The coexistence of the e-h and e-e orders is possible because 
e and h have opposite charges, whereas analogous phases are impossible in the case of two fermion 
species that have the same charge or are neutral. Most strikingly, the e-h order enhances the 
superconducting e-e order parameter by more than one order of magnitude as compared with that 
given by the BCS formula, for the same value of the effective e-e attractive potential X^"^ . This new 
phase should be observable in an e-h system created by photoexcitation in doped semiconductors 
at low temperatures. 



It is expected that electron- hole {e-h) systems created 
through photoexcitation of semiconductors exhibit var- 
ious phases depending on the material parameters and 
the densities iV^ and N*^ of electrons and holes, respec- 
tively ||l|] . Some of the interesting phases have been suc- 
cessfully observed, including the e-h liquid [gj and the 
Bose-Einstein condensation of excitons ||,Q| . These suc- 
cesses are largely due to the careful control of both the 
material parameters (by choosing a semiconductor) and 
the e-h density (through the excitation intensity). By 
further increasing the e-h density, one should observe an 
"e-/i BCS phase" at low temperatures, which is charac- 
terized by a nonzero e-h order parameter |q|. Although 
these phases are realized when TV^ = A^'*, one can also 
use doped semiconductors. In this case, one can control 
A^'^ and A^'' independently; N" - N'^ = N"^ through the 
donor (or acceptor) density A^**, and A^^ -I- A^'' through 
the excitation intensity ||5|]. The additional parameter 
A^^ may lead to a new quantum phase(s) at low temper- 
atures. When < |A^^| < A^" « A^'*, in particular, we 
may expect a new superconducting phase, which we call a 
multiply ordered superconducting (MS) phase, in which 
doped electrons form Cooper pairs moving in the e-h BCS 
state. On the other hand, most of the previous studies on 
superconductivity in the doped e-h BCS state (or related 
systems such as doped excitonic insulators and doped ex- 
citon systems) treated either the case where the doped 
electrons are located in a third band that is different from 
the bands of e-h pairs P,|7| , or the case A^'' <S A^*^ (or 
equivalently, A^*^ ^ A^'*)^Q. In these theories, the e-h 
BCS condensates (or excitons) work as polarizable me- 
dia, which induce a large attractive interaction between 
electrons, whereas the superconducting order parameter 
A"^*^ as well as the superconducting transition tempera- 
ture T^'^ is essentially given by substituting X'^'^ (effective 
dimensionless e-e attraction) and a cutoff parameter into 
the BCS formula or the McMillan formula |llO|. However, 
we expect the new MS phase when excess electrons (or 
holes) are doped in the e (/i) band where electrons (holes) 



forming the e-h BCS state are located., and when the e-h 
pairing is dominant, i.e., when < [A^^l <IC A^^ « A"'' 



and |Z\^ 



\A 



hh\ 



< \A 



eh\ 



Here, A'^'^ and A^'^ are the 



superconducting e-e and h-h order parameters, respec- 
tively, and A"^^ denotes the e-h order parameter. 

In this Letter, we explore the possibility of such a 
new MS phase, by studying the phase diagram of high- 
density e-h systems as a function of A^^. It is shown 
that the MS phase, which has the order parameters 
< {\A'"'l\A^^\) < |Z\'=''|, can be realized when < 
\N'''\ ^ N" K. N'\ Most strikingly, \A'"'\ is enhanced in 
the MS phase by more than one order of magnitude in 
comparison with the value given by the BCS or McMillan 
formula, for the same value of X^^ . 

We assume a three-dimensional, high-density, isotropic 
e-h gas at zero temperature, so that BCS-like mean field 
approximations work well 0. We decompose the inter- 
action Hamiltonian iJjnt into the short- and long-range 
parts, Hf^ and Hl^, respectively. The latter is related 
to long-range charge fluctuations, and will be discussed 
later. For the moment, we consider a charge-neutral 
region of unit volume, in which Hl"^ is irrelevant. In 
Hf^, the e-e, h-h and e-h interaction matrix elements 
are renormalized to effective values, U^"^ , U^^ and V^^ , 
respectively, due to the screening effect and negative (at- 
tractive) contributions from various bosonic excitations, 
such as lattice phonons and excitons. Since the bare 
value of V^^ is negative, the bosonic excitations make 
it more negative, whereas (basically positive) V^^ and 
U^^ are reduced. Hence, \U^^\ tends to be larger than 
|jjee|jjj/i/i|^ It was suggested that U"'' (and U^^) can be 
negative for some parameter values [^|j^,?J^. In the e-h 
BCS phase, there is an additional negative contribution 
from the Goldstone mode of the e-h order parameter. 
Since there is no reliable method of estimating V^^ and 
jjhh p^ j ^ ^g here treat them as given parameters, assum- 
ing that C/^^ U^^ < and \U''''\ > \U'"'\,\U''''\, and 
explore the phase diagram as a function of them and A^^ 
(> 0). We find that the minimal form of H^J^ for the MS 



state is 

fff„?= Y. '^kk'(4T'i-ki'i-k'iek'T+4i/iLkT^-k'Tek'i) 

k^tk' 

+ J2 ^kk'^ktetkie-kuek'T, (1) 

where eko- (/iko-) denotes the annihilation operator of e 
(h) with momentum k and spin a, and 



Utt = -V'\ 



kk' 

Tj-eh 
^kk' 



C/i 



kk' 



otherwise. 



(2) 



Here, V^'^, V^^ > are constants, and lUc (^ m) is a 
cutoff of the interactions Q . We assume that V^'^ > Y^^ 
so that \A^^\ :$> |Z\^'^|. Although we do not include an h- 
h interaction [/'*'' in H, we have confirmed that C/'''' can 
only modify the magnitude of the pair correlations by a 
factor of order unity, as long as \U'^'^\ < llf^^l. We do 
not include terms of the form eehh + h.c. (such terms are 
important in the case of two-band superconductors [Q ) , 
because in e-h systems intermediate processes involving 
such terms cost a large amount of energy of the order of 
the energy gap Eg (> {W'^l Hi). The total Hamiltonian 
without H^^ is denoted by H, and we put 
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FIG. 1. (a) n"(^k), (b) n''(Ck), (c) S^*", and (d) S^" plotted 
against ^k for solutions III (PU state; denoted by dotted lines) 
and V (MS state; solid lines), for TV ranging from 0.0 to 2.0 

0.25 and X"" = 0.10. 



in units of iV^pt. We take X^'' 



[(^k- t^)eLek^-F(Ck+ '^)hl^h^a] +Hsr. (3) 
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kCT 



Here, fi'' = Eg/2+n+i' and /i'' = Eg/2+n-i' [v > 0) are 
the chemical potentials of e and h, respectively, which are 
assumed to have the same energy dispersion k^/(2m) -|- 

J 



takes the same form as the BCS state, if (ek,o-, /j-k.-o-) 
is replaced by (ck,o-, C-k.-u)- This is the ordinary e-h 
BCS state |^ in nondoped {N'^ = TV'*) semiconductors. 
The energy cost of adding an electron-like quasiparticle 
to this state is £^k - i^, where Ei^ = y/£.'i + l^'^'^p- 
IH: |Z\'='*| 7^ 0, |Z\^"=| = 0; possible for smaU but finite 
N'^ > and v > |zi'^''|. Formally, this solution (whose 
wave function is denoted by |in)) is obtained from so- 
lution II (whose wave function |II)) by adding electron- 
like quasiparticles (whose annihilation operator eko-) up 
to E^ < u, i.e., |III) = (nL4.)|n), where JlL de- 



Eg/2. Moreover, N*^ _ z^k<T ^kcr^kCT, N^ _ z^ktr "-kCT^-kCT, notes the product over the range E^^ < v. Direct calcu- 

'nLeL'5'ka)|II)> where S'ko 



^k = k^/(2TO) — /i, and we take fi = 1 

We apply the mean field approximation that assumes 
the e-h correlation 5^ = (/i_k'xek't) a-^d the e-e corre- 
lation 5^ = (e-k'xek't)- We assume that (ft._kj,ek|) = 
(/i-kteki)- Using eq. @j, we find that the order pa- 
rameters, defined by Af = X^k' C^kk' ^\d and A'^ = 
Ek' t^kk' ^k' ) take simple forms: A<^ = A<'^ and A'^ = 



A^^ if lai < ujc 



^t 



otherwise. Here, A 



eh 



and A'^'^ are constants, which are taken to be real without 
loss of generality. We then diagonalize 7i, and obtain the 
self-consistent equations [|lj|. The system is character- 
ized by A''^ and dimensionless effective coupling constants 
X"^ = n-pV^^ , A*^^ = ufV^, where np is the density of 
states per spin at the Fermi surface. 

We have solved the self-consistent equations numeri- 
cally and found five solutions, which we denote by I-V: 
I: |Z\'='«| = |Z\ee| ^ 0; possible for all iV^. The one-particle 
distribution functions of e and h are n'^(^) = 0{i' — £,) and 
n''(.^) — 0{~v — ^), respectively, where 0{x) is the step 
function. 

II: \A''^\ ^ 0, |Z\'='=| = 0; possible for 7V^ = and j/ < 
|^eh|^ ^eh ^ for 1^1 < \A''% and the wave function 



lation shows that |III) 
annihilates an (ckm ^-k.-tr) pair. Therefore, n^{^) = 1, 
n^{C) = 0, and S^ = (e and h are unpaired) in the re- 
gion 1^1 < •\///2 — |Z\<='i|2 = ^^. This unpairing is demon- 
strated in Fig. n^ by dotted lines, which have discontinu- 
ities (secondary "Fermi surfaces") at ^ = ±^p. We call 
this solution the partially unpaired e-h BCS state (PU 
state). As N'^ increases, \A'^'^\ diminishes gradually until 
it vanishes at a certain value of A^^ , where this solution 
changes into solution I continuously. 
IV: |Z\'^''| = 0, iZi'^'^l / 0; possible for aU A^. Similar to 
solution I except that S^"" y^ ioi \^ - v\ < \A''''\. This is 
an ordinary superconductor of electrons. 
V: |Z\«''|, |Z\«'=| / QA^'^l » |Z\^«| because X"'' > X""); 
possible for small but finite A^^ > 0. Similar to solution 
IH except that S^^ ^ ii \^ ± Cpl < l^'^l, i.e., the e-e 
pair correlation exists around the secondary "Fermi sur- 
faces" (see solid lines in Fig. fh. This is the only solution 



where A*^ and A'^'^ coexist 1 15 . We call this solution the 



multiply ordered superconducting (MS) state. 

To identify what solution is physically realized. 
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FIG. 2. Plots oiE= (H) vs iV^ 
I-V, when X"'' = 0.25 and A^ 



N" - N'' for solutions 
0.20, for a fixed value of 
A'^ = (A'^'^ + A^ )/2. (The same plot is obtained for any values 
of N because of the normalization of both axes.) Aq'^ is the 
value of I /!'='' I at Af=" = of solution II. El and E^^ are the 
values of E at A'''^ — for solutions I and II, respectively. 
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FIG. 3. Phase diagram of a high-density e-h plasma for 
X"'' = 0.25. The MS phase is shadowed. lA""] is enhanced in 
the hatched region, where \A'"^\ takes maximum values with 
respect to N^ on the hne denoted as optimum doping. The 
vertical line at A'^ — corresponds to phase II. Dashed lines 
are guides for the eye (not calculated). 



compare their energies E = {H) for all values of N^ = 
N" - N'\ Note that we compare (H) rather than (H), 
because the natural parameters controlled directly by the 
photoexcitation intensity and the doping are N'^ and N^ 
rather than fj,"^ and //'' [g[. [If we compared {H), the 
discussion would be rather complicated because the re- 
lations between {N'',N^) and (^'^,/i'') are different for 
different solutions.] Figure shows i? as a function of 
N^, where N = (TV^ -I- 7V")/2. One recognizes that 
the solution V has the lowest energy for all values of 
N'^ /{ufAq^) < 1.86. However, care should be taken be- 
cause its curve is convex up, which might indicate a phase 
separation. We now show that the solution V is stable 
because e and h have opposite electrical charges. 

To show the stability, we first consider the unstable 
case where e and h denote some fermions that have 
charges of the same sign (including the neutral case). All 
calculations so far are also applicable to such a general 
case (as long as we regard U'^'^ and U'^'' as given param- 
eters). Since the energy of a single phase of V is larger 
than the average of the energies of two phases II and 
IV, the system should undergo a phase separation into 
two phases, one with the excess density N^ = < N^^^ 
(phase II) and the other with N^ > N^^^ (phase IV), 
where iV^^j denotes the average value of N^ of the total 
system. As N^^^ is decreased, the region(s) of phase IV 
becomes smaller, until the total system turns into a sin- 
gle phase II for N^^^ — 0. In a similar manner, one can 
also show the instability of phase III. This corresponds to 
the instability of the Sarma state |lq| that was discussed 
in the studies of superconductivity in a magnetic field. 

On the other hand, the situation is totally different if 
we return to the electron-hole system, where e and h have 
opposite charges —q and q {q > 0), respectively. If the 
phase separation occurred, each phase would have global 



net charge of density —q 6N = ~q{N'-^ — N"^) in phase II, 
and q5N in phase IV. The global charge would result in 
a large cost of the long-range part of the Coulomb energy 
_Elr = (Hf^), which we neglected in H, eq. (^). Taking 
i?LR into account, the total energy of such a nonuniform 
state would be Etot = J2i^i-^i + -^lR + -^b, where Vi 
and Ei denote the volume and the expectation value of 
H for phase i, respectively, and Eb the boundary en- 
ergy. Since Eb > 0, E^.t > J2^v,E, + ^lr = ^^ot- 
When the system is in a single phase V, then _Elr = 
and E'tot equals E of phase V of Fig. g. On the other 
hand, when the system is separated into cells of phases 
II and IV, one obtains a finite -Elr, which would be dom- 
inated by the electrostatic energy. We find that E^^^ is 
smaller than E^]^^^^^ (< EI^^^^), hence a single phase 
V should be realized, except for very small values of 

jyx jf ^eh ^ ]^^ jjg^g^ ^eh ^ „p/(7r|Z\e'«|) is the Pip- 

pard length of e-h pairs, which gives the minimum allow- 
able size of cells of phase II, whereas Lc ~ \/n/{q'^n-p) 
is the maximum allowable size of cells of phase II, be- 
low which the energy decrease by the phase separation 
overcomes the increase of -Elr, where k denotes the di- 
electric constant of the semiconductor. The condition 
L^^ 3> Lc is always satisfied at high densities, i.e., when 
r, EE [3/{ATTN)]^^^mq^/{ATrKh^) < 1 and N" « N'^ 
(w N). In fact, in this case we may assume a screened 
Coulomb interaction for V'^'^ and can easily show that 
L'^^/Lc > 10^ for any Tg ^ 1. Therefore, at high den- 
sities (rg ^ 1) the phase V is always stable against 
phase separation. Moreover, considering the large value 
10^, we may extend this conclusion to densities where 
Ts « 1 Jl^. The condition Ts < 1 can be realized when, 



e.g.. 



0.1m: 



0' 



IOkq and N = 10^ 



for which 



we obtain rg = 0.54, where TOq is the free electron mass 
and kq the dielectric constant of vacuum. 

By performing calculations of the energies of various 
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FIG. 4. The 7V^ dependence of \A''^\ in the MS phase for 
°'' = 0.25. The inset shows the enhancement factor Q vs 



solutions as a function of A*^' 
phase diagram shown in Fig. 



and A^^, we obtain the 
I Phase V (MS phase), 



for which \A \ ^ |Z\'^^| > 0, is reahzed in the shadowed 
region. This phase changes continuously into phase II as 
N^ -^ 0. On the other hand, this phase changes into 
phase IV as N'^ increases, because then the Fermi en- 
ergies of e and h become separate, which prevents e-h 
pairing. To study the N^ dependence in more detail, we 
plot lA"^] in the MS phase in Fig. |[ One sees that |Z\'='=| 
strongly depends on N'^ and takes a maximum value Ziopt 
at an optimum N"" (= iV^pJ, although N" = N + iV^/2 
is almost constant (because N^ ^ N) . This is in marked 
contrast with the single-component case, where the BCS 
theory gives \A^^\ = 2uJc exp{—l / X^^) = Aj^^g, which 
is constant when N^ and X^^ are constant. To com- 
pare /iopt with Z\g'^(^g, we plot the enhancement factor 
= Z\^^j/Z\|'cs in the inset of Fig. |. Since zi|^cs is 
the magnitude of \A'^'^\ which would arise if there were 
no e-h interaction, the fact Q > 1 implies that \A'^'^\ 
is enhanced by the presence of \A'^'''\. It is found that 
an enhancement by one order of magnitude or larger oc- 
curs as A"^^ becomes smaller. As X^^ -^ 0, in particular, 
[A'^'^l diminishes much more slowly than A^q^, resulting 
in a very large enhancement factor Q. Since the super- 
conducting transition temperature T^"^ oc |/i'^'^|, the MS 
phase should have a T^"^ which is much higher than the 
BCS transition temperature p8| . The parameter region 
in which this enhancement occurs is hatched in Fig. 0, 
where the optimum doping line is also shown. 

To reveal the physical mechanism leading to the en- 
hancement of |/i'^'^|, we reexamine the e distributions of 
the MS and PU states of Fig. |l](a), which is schematically 
plotted in Fig. |. In the PU state, for which lA""] = 0, 
excess electrons are concentrated in the region |^| < fp, 
and the secondary (lower and upper) "Fermi surfaces" 
appear at ^ = ±^p. In the MS state, electrons form e-e 
pairs to benefit from the attractive e-e interaction and 
consequently S^ develops around those "Fermi surfaces" 
[Fig. |^(d)]. However, this pairing necessarily broadens 
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FIG. 5. Schematic illustration of the enhancement mecha- 
nism. The electron distributions n'^(^) for the PU state (dot- 
ted line) and the MS state (solid line) are shown. 



the "Fermi surfaces" , resulting in increase of the kinetic 
energy E'kin- Hence, the total energy is minimized for 
an optimum S^. Although this mechanism is similar to 
the BCS instability in ordinary superconductors, there 
is an essential difference: In the BCS state, the broad- 
ening of the Fermi surface of width 6^ increases the ki- 
netic energy by SE^-in ~ np|(5^p « np|Z\'='=p. In the MS 
state, on the other hand, the broadenings of the lower 
and upper "Fermi surfaces" change the kinetic energy by 
SE<^ « -nplA'^'l^ and SE>^ « +nF\A'=''\^, respectively 
[when X'^^ (thus \S£^\) is small]. These two contributions 
cancel to give a relatively small increase in the kinetic en- 
ergy (JE'kin = i^^kin + '^-^kin- ^his is the Origin of the huge 
enhancement factor Q for small A*^"^. Since the cancella- 
tion is not perfect, \A'^'^\ of course takes a finite value. 
As X^^ is increased, the cancellation becomes less com- 
plete, and Q is reduced. The existence of an optimum 
doping can be understood in a similar manner: If N^ 
is too small, no room for broadening is left between the 
two "Fermi surfaces" . Conversely, if A^^ is too large, the 
magnitude of the jump of n'^{^) at the two "Fermi sur- 
faces" differs, which makes the cancellation of contribu- 
tions from the two "Fermi surfaces" less perfect, result- 
ing in a smaller jZi'^'^l. Therefore, \A'^'^\ takes a maximum 
value at an intermediate value of N'^ . 

Finally, we point out that although we have assumed 
the e-h BCS state in semiconductors, the same mech- 
anism for the enhancement of the e-e correlation upon 
doping may be expected also in other ordered phases of 
other materials, if the electron distribution without the 
e-e correlation is similar to the dotted hue of Fig. g. 
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